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Abstract. We prove that for a finite collection of real- valued functions fi, . . . , f„ on the group 
of complex numbers of modulus 1 which are derivable with Lipschitz continuous derivative, the 
distribution of (tr/i, . . . ,trf„) under the properly scaled heat kernel measure at a given time on the 
unitary group U(A'') has Gaussian fluctuations as A'^ tends to infinity, with a covariance for which 
we give a formula and which is of order A'^^^. In the limit where the time tends to infinity, we 
prove that this covariance converges to that obtained by P. Diaconis and S. Evans in a previous 
work on uniformly distributed unitary matrices. Finally, we discuss some combinatorial aspects of 
our results. 
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1. Introduction 

In [7], P. Diaconis and S. Evans studied the fluctuations of the trace of functions of a unitary 
matrix picked uniformly at random. Let us recall briefly their main result. If [/ is a unitary matrix 
of size > 1 and / a real- valued function on the set U of complex numbers of modulus 1, then the 
eigenvalues Ai, . . . , Aat of U belong to U and tr/(C/) = jj^ J2iLi f{^i)i where tr is the normalized 
trace (so that tr(/jv) = 1) and the matrix f(U) is obtained from U and / by functional calculus. 
By rotational invariance of the Haar measure, it is easy to see that if / : U — > M is defined almost 
everywhere, integrable and has zero mean on U then tr f(U) is defined for almost every U and, seen 
as a random variable under the Haar measure, also has zero mean. 

The function / being fixed, tr/ can be seen as a random variable on the unitary group U(A^), 
endowed with the Haar measure, for all > 1. Thus, the single function / gives rise to a sequence 
of random variables indexed by the integer N, which is their main object of study. In order to 
understand the behavior of this sequence, a fundamental fact, which has been proved and used 
extensively in this context in [7j, is the following: for all p,q G Z, one has E[tr(C/^)tr(C/'J)] = 
6p^qN~'^ mm{p, N). Using this, one can easily check that, if / is square-integrable on U, then the 
variance of tr/ converges to as tends to infinity. Moreover, if / belongs to the Sobolev space 
H2(iJ) (see Definition 18.11 below) . then the series of the variances of tr/ on U(A^) converges, which 
gives a strong law of large numbers. 

The main result of [7] is that the fluctuations of tr/ under the Haar measure are asymptotically 
Gaussian. More precisely, they have proved that if / belongs to i?2(U) and has zero mean on U, 
then Ntrf converges in distribution to a centered Gaussian random variable with variance equal 
to the square of the ffz-norm of / (see Theorem 18.21 below for a precise statement). 

In this paper, we consider the fluctuations of tr/ when the unitary matrix is picked not under the 
Haar measure, but rather under the heat kernel measure at a certain time. The heat kernel measure 
at time T is the distribution of Un{T), where (f/Ar(t))t>o is the Brownian motion on U(A^) issued 
from the identity matrix, that is, the Markov process whose generator is the Laplace-Beltrami 
operator associated to a certain Riemannian metric on U(A^). The choice of a Riemannian metric 
that we make is explicited at the beginning of Section [2j Apart from being one of the most natural 
stochastic processes with values in the unitary group, the Brownian motion arises for example in 
the context of two-dimensional U(A^) Yang-Mills theory (p ^ [TO l [9]). 

Let / : U ^ M be a function, as above. Once a time T > is fixed, tr/ is a random variable 
on U(A^) for each A'^ > 1, the unitary group being endowed with the heat kernel measure at time 
T. With our choice of Riemannian metric, it is known since the work of P. Biane [2] that if / is 
continuous, then tr/ converges almost surely towards the integral of / against a probability measure 
on U, which is characterized by the formula (jH) below. By this almost sure convergence, we 
mean that the expectations of these variables and the series of their variances converge. For all 
T > 0, the measure vt is absolutely continuous with respect to the uniform measure on U, with 
a density which unfortunately cannot be expressed in terms of usual functions. Its support is the 
full circle only for T > 4. For T £ (0,4), its support is an arc of circle containing 1, symmetric 
with respect to the horizontal axis, which grows continuously with T, and for the width of which 
a simple explicit formula exists. In fact, as A^ tends to infinity, not only the distribution of the 
eigenvalues oi (T) but the Brownian motion itself as a stochastic process converges in a certain 
sense towards a limiting object called the free multiplicative Brownian motion, which is defined in 
the language of free probability. The measure i't is the non-commutative distribution of this free 
process at time T and can be considered as a multiplicative analogue of the Wigner semi-circle law. 
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The main result of this paper is that for any function / : U ^ M with Lipschitz continuous 
derivative, the fluctuations of iVtr/ are asymptotically Gaussian with variance /), where 

ax is the quadratic form defined in Definition 12.41 This definition of crx{f,f) involves three free 
multiplicative Brownian motions which are mutually free and the functional calculus associated to 
/'. It makes sense for functions of class C^, or at best for absolutely continuous functions. An 
alternative definition of o"t(/, /) is given by Definition 18. 101 in terms of the Fourier coefficients of / 
and the solution of an infinite triangular differential system (see Lemma [8.7p . We prove that, when 
T is large enough, this second definition makes sense for functions in the Sobolev space H^ll]")^ 
which are not even necessarily continuous. 

Moreover, we prove that, as T tends to infinity, crT{f,f) converges towards the square of the 

2-norm of /. This convergence is consistent, at a heuristic level, with the result of P. Diaconis 
and S. Evans, since the Haar measure is the invariant measure of the Brownian motion, and its 
limiting distribution as time tends to infinity. 

For small values of T, the analysis seems much harder to perform. We have no expression of the 
covariance other than Definition 12.41 and it seems plausible, considering the limiting support of the 
distribution of the eigenvalues of C/at (T) and some puzzling numerical simulations (see Figure [T] in 
Section [8]) , that the largest space of functions / for which Ntrf has Gaussian fluctuations might 
depend on T, say for T < 4. Unfortunately, we have no precise conjecture to offer in this respect. 

The understanding of global fluctuations of random matrices has been widely developed in the 
literature using various techniques. By combinatorial methods applied to the computation of 
moments, Ya. Sinai and A. Soshnikov [25] derived a central limit theorem (CLT) for moments of 
Wigner matrices growing as o{N^^^). An important breakthrough is the work of K. Johansson |15j 
where he got, using techniques of orthogonal polynomials on the explicit joint density of eigenvalues, 
a CLT for Hermitian or real symmetric matrices whose entries have joint density e^*''^^*^), for a 
large class of potentials V. Recently, M. Shcherbina [24] has been able to lower, in the symmetric 
case, the regularity of those functions for which the CLT holds. The study of Stieltjes transform 
for this purpose, initiated by L. Pastur and others |221l23j. has recently given some striking results, 
among which one can cite the works of G. W. Anderson and O. Zeitouni [T] or W. Hachem, P. 
Loubaton and J. Najim [13]. Recently S. Chaterjee [5] proposed "a soft approach" based on second 
order Poincare inequalities. 

Some tools of free probability will play a key role in our analysis. In this framework, a notion 
of second order freeness was developed in [211 ISOl [6]. In particular, the second paper [20] of the 
series deals with unitary matrices and the results that it contains might be relevant to the problem 
under consideration. 

The technique of proof that we have chosen is rather of the flavour of the one introduced in 
[3]. Therein, T. Cabanal-Duvillard proposed an approach based on matricial stochastic calculus 
to get a CLT for Hermitian and Wishart Brownian motions but also for several Gaussian Wigner 
matrices. In this direction we can also mention a CLT for band matrices obtained by A. Guionnet 

The paper is organized as follows : Section [2] is devoted to defining the Brownian motion on the 
unitary group, recall from [2] its asymptotics, define the proper covariance functional and state our 
main result (Theorem 12. 6p . In Section [3l we present the structure of the proof of our main theorem 
by introducing a family of martingales (see Equation Q) that will be the main object of study. 
The proof will in fact boil down to proving the convergence of the bracket of these martingales 
(Section [5]) and to controlling the variance of this bracket (Section [6]), relying on some technical 
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results on the functional calculus on U(A^) gathered in Section [H In Section [71 we extend our result 
to other Brownian motions on the unitary group and to the Brownian motion on the special unitary 
group. Section [5] is devoted to the study of the covariance for large time, in connexion with the 
CLT for Haar unitaries [7]. Finally, in Section [U we discuss a combinatorial approach to some of 
our previous results. 



2. The Brownian motion on the unitary group 

2.1. The stochastic differential equation. Let > 1 be an integer. We denote by U(A'^) the 
group of unitary N x N matrices and by u(A^) its Lie algebra, which is the space of anti-Hermitian 
N X N matrices. We denote by 1^ the identity matrix. We will use systematically the following 
convention for traces: we denote the usual trace by Tr and the normalized trace by tr, so that 
Tr(/Ar) = N and tr(/Ar) = L 

Let us endow u(iV) with the real scalar product (X,y)u(jv) = NTy{X*Y) = -NTy{XY). We 
denote by || • ||u(Ar) the corresponding norm. 

The scalar product ( • , • ) u{n) determines a Brownian motion with values in u( A^) , namely the 
unique Gaussian process {Kiy{t))t>o with values in u(A^) such that 

Vs, t>0,yA,B£ n{N), E[{A, i^v(s))u(JV) {B, A',v(t))u(7V)] = min(s, t){A, S)^^). 

Equivalently, let C^i, -Dfc)A;,i>i be independent standard real Brownian motions. Then K]\f(t) 
has the same distribution as the anti-Hermitian matrix whose upper-diagonal coefficients are the 
'^^{Bki{t) + iCki{t)) and whose diagonal coefficients are the -^Dk{t). 

The linear stochastic differential equation 

dUNit) = UN{t)dKN{t) - \uN{t)dt (1) 

admits a strong solution which is a process with values in M7v(C). This process satisfies the identity 
d{UNU'^){t) = 0, as one can check by using Ito's formula. Hence, this equation defines a Markov 
process on the unitary group U(A^), which we call the unitary Brownian motion. The generator 
of this Markov process can be described as follows. Let (Xi, . . . ,^^2) be an orthonormal basis 
of u(A^). Each element X of u(A^) can be identified with the left-invariant first-order differential 
operator Cx on U(A^) by setting, for all differentiable function F : 1}{N) M and all U E U(A^), 

iCxF)iU) = ^ F(C/e*^). (2) 
dt\t=o 

The generator of the unitary Brownian motion is the second-order differential operator 

2 2^ 

k=l 

This operator does not depend on the choice of the orthonormal basis of u(A^). We denote the 
associated semi-group by {Pt)t>o- From now on, we will always consider the Brownian motion 
issued from the identity matrix, so that Un{0) = In. 

The stochastic differential equation satisfied hy Un can be translated into an Ito formula, as 
follows. 
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Proposition 2.1. Let F : R X \J{N) ^ R be a function of class C^. Then for all t > 0, 

F{t,UN{t)) = F{0,lN) + y] [ {Cx,F){s,Un{s)) d{Xk,KN)^iN){s) 

Jo 



k=l ■ 

+ (^^AF + dtF^{s,UN{s))ds, (3) 

and the processes {(X^, i^7v)u(Af) • ^ ^ {!)••• j-^^}} ^'^c independent standard real Brownian mo- 
tions. 

This result is classical in the framework of stochastic analysis on manifolds (see for example 
|14j). but since our whole analysis relies on this formula and for the convenience of the reader, we 
offer a sketch of proof in this particular setting. 

Proof. For all a,b £ {1, . . . , N}, let £ab '■ MAr(C) C denote the coordinate mapping which to a 
matrix M associates the entry Mab- Let also dab denote the partial derivation with respect to the 
a6-entry. The definition of Cx given by ([2]) makes sense for any matrix X. One can check the 
following identities: 

N N 

MX G MAr(C), Cx= SacXcbdab and Cx - Cx^ = ^ SacXcbSa'c'Xc'b'dabda'b' , 

a,b,c=l a,b,c,a' ,b' ,c'=l 

7V2 N 

A = ^C + '^ ^ £ac{Xk)cb£a'c'iXk)c'b'dabda'b'i 
k=l a,b,c,a',b',c'=l 

where C = Ei=i^f- Mor eover, C = —In, regardless of the choice of the orthonormal basis 

{Xl, . . . , X^2). 

Any smooth function F -.Rx U(A^) is the restriction of a smooth function defined on M x MAr(C). 
Applying the usual Ito formula to this extended function and using the identities above leads im- 
mediately to ([3]). □ 



2.2. The free multiplicative Brownian motion. We are interested in the large behaviour 
of the stochastic process Un issued from 1^. Philippe Biane has described in [2] the limiting 
distribution of this process seen as a collection of elements of the non-commutative probability 
space (L°° ® MAr(C), E (8) tr). We start by describing the limiting object. As a general reference on 
non-commutative probability and freeness, we recommend |26j . 

Definition 2.2. Let {A,t) be a (non- commutative) * -probability space. A collection of unitaries 
iut)t>o in -4 is called a free multiplicative Brownian motion if the following properties hold. 

1. For all < ti < . . . < tn, the elements ut^ , lifjU^^ , . . . , ut„Uf^_^ are free. 

2. For all < s < t, the element utu* has the same distribution as ut-s- 

3. For all t > 0, the distribution of ut is the probability measure ut on \} = {z £ C : \z\ = 1} 
characterized by the identity 

/ ^ . dvt{i) = l + z, (4) 

1 _ ^e*^e^e 

valid for z in a neighbourhood of 0. 

The following result has been proved by P. Biane. The second assertion follows from the first by 
a general result of D. Voiculescu. 
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Theorem 2.3. The collection {U]\f{t))t>o of non-commutative random variables converges in dis- 
tribution, as N tends to +00, towards a free multiplicative Brownian motion. 

Moreover, if U^^ , , ■ ■ ■ , are n independent sequences of unitary Brownian motions, the 
family {{U^\t))t>o, {^n^ {'t))t>o, ■ ■ ■ , {UN^\'t))t>o) converges in non- commutative distribution, as 
N tends to infinity, towards {{Uf^^)t>o, {u[^'^)t>o, . . . , iui^^)t>o) where u''^\ . . . jU^*^) are n free mul- 
tiplicative Brownian motions which are mutually free. 

2.3. Statement of the Central Limit Theorem. Recall that U denotes the group of complex 
numbers of modulus 1. Let / : U — > M be a function. Then, by the functional calculus, / induces a 
function, still denoted by /, from U(A^) to M7v(C). Moreover, for all unitary matrix U , the matrix 
f{U) is Hermitian. 

We endow U with the usual length distance, that is, the distance such that d(e*",e*^) = \a — I3\ 
for all a,/3 G M such that |q; — /3| < vr. Accordingly, we define the Lipschitz norm of a function 
/ : U ^ R as follows: 

|/(z)-/H| 

||/||Lip = sup ^ . 

z,w&j,ZTtw a[z,w) 

Note that, if / is Lipschitz continuous and z, w belong to U, then the following inequalities hold: 
\f{z)-f{w)\ < \\f\\ur,d{z,w) < |||/||Lip|z-H. 

By the derivative of a differentiable function / : U ^ M, we mean the function /' : U — > M 
defined by 

V. . U, /'(.) = lim IM^IZIMI. 

' ' ^ ' h~*Q h 

We denote by L^(U) the space of integrable functions on U, with respect to the Lebesgue measure. 
We denote by C'^{V) the space of continuously differentiable functions and by C^'^(U) the subspace 
of C^(U) consisting of those functions whose derivative is Lipschitz continuous. We define a family 
of bilinear forms on C^iJ]) as follows. 

Definition 2.4. Let [A, r) be a C* -probability space which carries three free multiplicative Brownian 
motions u,v,w which are mutually free. Let T > be a real number. Let /,(/ : U ^ M &e two 
functions 0/ C^(U). For all s G [0,T], we set aT,s{f,g) = T{f'{usVT-s)g'{usWT-s))- Then, we 
define 

(^T{f,g)= CrT,s{f,g)ds= T{f'{UsVT-s)g'{UsWT-s)) ds. 



Lemma 2.5. For all T >0, ax is a symmetric non-negative bilinear form on C^(U). 

Proof. The symmetry of ax comes from the fact that the triples {u, v, w) and {u, w, v) have the 
same distribution. In order to prove the non-negativity, let us realize {u, v, w) on the free product 
of three non-commutative probability spaces. So, let {Au,Tu), {Av,Ty) and {Aw,Tw) be three non- 
commutative probability spaces which carry respectively u, v and w. We consider their free product, 
so we define A = Au * Av * Aw and t = Tu * Ty * Tw We also use the notation Tu, Ti,,Tw for the 
partial traces on A. Then 

(^T{f,f)= / Tu{Ty{f'{UsVT-s))rw{f'{UsWT-s))) ds = / Tu{Ty{f' {UsVT-s)f) ds > 0, 







the positivity coming from the fact that f {ugVT-s) is self-adjoint. □ 



We will use the notation o"r(/) = o"t(/, /)• Let us state our main result. 
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Theorem 2.6. Let T > be a real number. Let n > 1 be an integer. Let /i, . . . , : U M 

be n functions of C^'"^(U). Let us define a n x n real non-negative symmetric matrix by setting 
^r(/i) • • • ) /n) = (<7r(/i) /j))jje{i,...,n}- Then, as N tends to infinity, the following convergence of 
random vectors in M" holds in distribution: 

N{tTf,{UN{T))-E[tTf,{UN{T))]\^^,_^y AA(0,St(/i,...,/„)). (5) 



3. Structure of the proof 

For r = 0, the result is straightforward. Let us choose once for all a real T > 0. In order to 
study the left-hand side of (|14p . we write each component of this random vector as the difference 
between the final and the initial value of a martingale. To do this, let {J^N,t)t>o denote the filtration 
generated by the unitary Brownian motion C/jy- To each function / of L^(U) we associate a real- 
valued martingale (MjY(t))(g[o^T] by setting 

M{,{t)=E[tTf{UN{T))\J^N,t]. (6) 
The left-hand side of (jl4p is simply N (m(i{T) — Mh{0)] and we are going to study the 

V /ie{l,...,n} 
f . 

quadratic variations and covariations of the martingales Mjy\ In order to state the main technical 
results, let us introduce some notation. 

Recall that the gradient of a differentiable function F : U(A^) ^ C is the vector field on U(A^) 
defined by VF = J2k=i(^XkT')Xk, where {Xi, . . . ,Xj^2) is an orthonormal basis of u(A^). To each 
pair of functions f,g& -^^(U) we associate a function sj^^ on [0,r) x U(A^) by setting 

eI^'^{s,U) = iV2(V(Pr-.(tr/))(C/),V(PT-.(tr5))(t/))u(^). 

Let us check that this function is well-defined. By the Weyl integration formula, the fact that / 
is integrable on U implies that tr/ is an integrable function on U(A^). Hence, for all s G [0,T), 
Pr-s(tr/) is a function of class C°° on U(A^) and eJ^^ is well defined. 

Proposition 3.1. Consider f,g^ L^(U). With the notation introduced above, the following prop- 
erties hold. 

1. For all t G [0,T], the quadratic covariation of the martingales NAL^ and NMf^ is given by 

{NM^^,NMf,)t= f E^j^3{s,UN{s))ds. 
Jo 

2. Assume that f and g are Lipschitz continuous. Then for all s G [0, T) and all U £ U(A^), 
\Ej^^ {s,U)\ < (||/||Lip+ IbllLip)^- Moreover, if f and g belong to C^(U), then the following 
convergence holds: 

3. Assume that f and g belong to C^'^{U). Then the following estimate holds: 

sup Var(F^f(s,C/^(s))) = 0{N-^). 
se[o,T) 

Let us show that these results imply Theorem 
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Proof of Theorem \2.6l For all > 1, define a M"-valued martingale Qn = (Qjv • • • ' Q%) setting 

It is a martingale indexed by [0,r], issued from and 



QN{t) = N(Ml^{t)-M^^{0)) 



with the same bracket as N ( M-j^ 



je{i,...,n} 



For alU = (Ci, . . . , Cn) G IK" and all t £ [0, T], set 



/ n ^ n t \ 

RN{t) = exp iy]^jQ]v(t) + o S.jCk / o-T,s(/j, /fc) (is • 



Ito's formula yields 



j,k=i 



Thus, 



||t||2 ^iT^ II W| 

\E[RN{t)-l]\ < ' . = l...nl'^^' 



2 J 
For fixed j and A;, the last integral is smaller than 



max E /* aT,s{fj,fk)-E^^'^\s,UN{s)) 

,k=l...n Jq 



ds. 



aTAfjJk)-E\Ey''{s,UN){s) 



ds + E 



Ej^'^''{s,UN){s)-E 



E^^'^''{s,UN)is) 



ds. 



By the second part of Proposition 13.11 and by the dominated convergence theorem, the first 
integral tends to as tends to infinity. The square of the second integral is smaller than 
t Jq Var(£'^' Ui\i{s))) ds, which, thanks to the third part of Proposition 13. II and by dominated 
convergence again, tends also to 0. Finally, we have proved that 



G M", lim E 



which, for t = T, yields the expected result 



1 



j,k=l 



□ 



In Section[31 we collect some technical results that we use in Sections [5] and[6]to prove Proposition 

EH 

4. Regularity of the functional calculus 

In this section, we relate the regularity of a function / : U ^ M to the regularity of the functional 
calculus mapping / : U(A^) M7v(C) and the function tr/ : U(iV) M. We start with a result 
which, logically speaking, is not necessary for our exposition, but which is the simplest instance of 
a crucial phenomenon. 

4.1. Lipschitz norms. The group U(A^) becomes a metric space when it is endowed with the 
Riemannian distance, denoted by d, associated to the Riemannian metric induced by the scalar 
product (•, •)u(7v) on u(A^). We denote by ||-F||Lip the corresponding Lipschitz norm of a function 
F : U(iV) ^ M, that is, 

Ili^llLip = sup I ■■U,V£ U(iV), [/ / f} . 
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Proposition 4.1. Let / : U ^ M 6e a Lipschitz continuous function. Then trf : U(A^) — > M is also 
Lipschitz continuous and 

||tr/||Lip = ^||/||Lip. 

Note that this result can be compared to Lemma 1.2 in [12], where it was a key point towards the 
concentration resuhs for Wigner and Wishart random matrices. In order to prove this proposition, 
we use the fohowing lemma. 

Lemma 4.2. Let U and V be two elements of\J{N). Then there exists A,B£ U(A^) such that 
AUA-^ and BVB~^ are diagonal and d{AU A'^ , BV B'^) < d{U,V). 

Proof. Let O be the conjugacy class of V. It is a compact submanifold of U(A^). Let V' be a point 
of O which minimizes the distance to U. Let 7 : [0, 1] — > U(A^) be a minimizing geodesic path from 
V to U parametrized at constant speed. It is thus of the form ^(t) = V'e^^ for some Z G u(A^). 
Since V' minimizes the distance to U, the vector 7(0) is necessarily orthogonal to the tangent space 
TyO. This space TyO, identified with a subspace of u(A^) by a left translation, is the range of the 
linear mapping Ad{V'~^) — Id. Hence, Z belongs to the kernel of the adjoint linear mapping, that 
is, to the kernel of Ad{V') — Id. In other words, V'ZV'~^ = Z. It follows that Z and V can be 
simultaneously diagonalized, in an orthonormal basis, and the same is true for V' and V'e^ = U. 
Finally, V' and U are conjugated by a same unitary matrix to two diagonal unitary matrices. The 
result follows easily. □ 



Proof of Proposition \4-.l\ Let / : U ^ M be Lipschitz continuous. Consider U and V in U(A^). 
Thanks to Lemma 14.21 let us choose U' and V' which are both diagonal, conjugated respectively 
to U and V, and such that d{U' ,V') < d{U,V). Let us write U' = diag(e*"i , . . . , e*"^) and 
V = diag(e*^i, . . . , e*^^) in such a way that \(3j — aj\ < vr for all j G {1, . . . , N}. Let us compute 
d{U', V). It is equal to d{lN, U'~^V'), hence to 



d(/;v,e''''''^^^^-"^'-'^^-"^)) = ||idiag(/?i - ai, . . . ,/3^ - aN)L(N) 



N 

2 



It follows that d{U, V) > X^^i — cxjl- On the other hand, 

TV N 

|tr/(y)-tr/(C/)| <-j;|/(e^/^0-/(e^"OI <^ll/llLipEl/^i-«^-| ^]v''^ 

j=i j=i 

This proves the inequality ||tr/||Lip < ;^||/||Lip- By choosing a, (3 such that |/(e*^) — /(e*")| is close 
to ||/||Lip|/? — a| and by considering U = c^^/at, V = e^^L^, one verifies that the opposite inequality 
holds. □ 



Let us make a short heuristic comment on this result. The scalar product which we have chosen on 
u(A^) corresponds to a metric structure on U(A^) which gives this group the diameter d{L]y, —/at) = 
||idiag(7r, . . . , vr)||u(7v) = ^'^i of the order of N. The function / : U — > M being fixed, the variations 
of the function trf : D{N) — > M are of the same order of magnitude as those of / but occur on a 
space N times as large. This makes the equality that we have juste proved plausible. 

In the same order of ideas, note that the distance to the origin at time T of a linear Brownian 
motion in a Euclidean space of large dimension d is, by the law of large numbers, of the order of 
\/ dT. Assuming that the Brownian motion on the unitary group behaves in a comparable way, and 
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considering the fact that the dimension of U(A^) is N'^, this indicates that the Brownian motion 
Un{T) might be at a distance of order of 1^, thus a fraction of the diameter of U(A^) which 

does not depend on N. This gives an intuitive justification for the choice of the normahzation. 

4.2. First derivatives. We are now going to prove that the functional calculus induced by / is 
differentiable when / is differentiable, and to compute its differential. For this, we introduce some 
notation. Let / : U ^ C be a differentiable function. Let us define a function D/ : U x U ^ C by 
setting 

C f{z)-f{w) ;f ^ 

\fz,w€V,Dfiz,w) = { f-f !'/^^' 

The function D/ is symmetric and, if / is C^{U), it is continuous and bounded by ^||/'||oo- Note 
that D/ takes its values in C even if / is real-valued. 

If the function / is only Lipschitz continuous, then it is differentiable with bounded differential 
outside a negligible subset of U, and the definition of D/ still makes sense outside the corresponding 
negligible subset of the diagonal of U x U. Moreover, outside this subset, the inequality | D/(z, w) \ < 
f ll/'lloo holds. 

If i7 is a unitary matrix, we denote by Ljj and Ru the linear operators on MAr(C) of left and 
right multiplication by U respectively. These operators commute and they are normal with respect 
to the scalar product {A, B) = NTh:{A* B) on M7v(C). In fact, = L^-i and = Ru-i. Hence, 
if (7 is a function on U x U, then g{Lif, Rjj) is a well-defined endomorphism of M7v(C). Even when 
/ is only Lipschitz continuous, Df{Lu,Rij) is well-defined for almost all U £ U(A^). 

Let us define a special orthonormal basis of u{N). We use the notation (-EjA:)j,fce{i,...,7V} 
the canonical basis of Mjv(C). For all j,k with I < j < k < N, set Xjk = {Ejk — Ej^j) 
and Yjf; = --^=[Ejj. -\- Ej^j). For all j G {1, . . . , A^}, set Hj = -j=Ejj. These matrices form an 
orthonormal basis of u(A^). 

Proposition 4.3. Let f : V ^ C be a differentiable function. Let U be an element of U(iV). Let 

X be an element ofu{N). Then 

^ f{Ue''') = {Df{Lu,Ru)){UX). (7) 
dt\t=o 

In particular, when U is a diagonal matrix with diagonal coefficients {ui, . . . ,un), the following 
equalities hold. 

1. For all 3 G {1, . . .,N}, f^^^^J (C/e*^^) = Df{uj,uj)UHj. 

2. For all j,ke {!,..., N} with j < k, f^^^^^f {Ue'^^^) = Df{uj, Uk)UXjk. 

3. Forallj,k€{l,..., N} with j < k, f^^^^J (t/e*^^'=) = Df{u,,Uk)UYjk. 

If f is only Lipschitz continuous, then the same conclusions hold for almost all U £ U(A^). 

Proof. We will give the proof under the assumption that / is differentiable. The extension to the 
Lipschitz continuous case is straightforward. Let us start by proving the part of the statement 
which concerns a diagonal matrix U. 

1. Since Ue*^^ is diagonal, this assertion is proved by an easy direct computation. 

2. This case is less trivial. Let us assume that Uj ^ u^. Then for small t, there is a unique 
pair of continuous functions {uj{t),Uk{t)) such that the spectrum of Ue*-^^'' is deduced from that 
of U by replacing Uj and Uk respectively by Uj{t) and Uk{t). The functions Uj and ut are in fact 
smooth and they satisfy u'j{^) = ^^^(0) = 0, an equality which can be phrased by saying that the 
right multiplication by e*^^* does not affect the spectrum of U at the first order. 
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Let D{t) be the diagonal matrix obtained from U by replacing Uj and Uk by Uj{t) and Uk{t) 
respectively. By diagonalizing Ue''^^'' for small t, one can find a unitary matrix P{t) which depends 
smoothly on t, such that P{0) = In, such that the only non-zero off-diagonal terms of P{t) are 
P{t)jk and P{t)kj, and finally such that 

Ue^^^" = P{t)D{t)P{t)-^. (8) 

By differentiating with respect to t at t = 0, one finds 

UXjk = [P'{0),U], 

from which one deduces that P'{0)jk = ^ and P'{Q)kj = ^^X^ ■ By applying / to 

both sides of ([H) and then differentiating again with respect to t at t = 0, we find 

Knowing the off-diagonal terms of P'{0) is enough to compute this bracket and we find the expected 
result. The case where Uj = is left to the reader, as well as the third assertion. 

Let us now turn to the first part of the statement, where no assumption is made on U . Let us 
first prove that ([7]) is true when C7 is a diagonal matrix with diagonal coefficients (ui, . . . , un). 

In this case, for all j, E {!,..., A^}, the matrix Ej^ is an eigenvector for Ljj and Ru, with the 
eigenvalues Uj and respectively. Hence, by definition of D/, Ej^ is an eigenvector of Df{Lu, Ru) 
with the eigenvalue Df{uj,Uk). The validity of ([7]) in this case follows, because UHj (resp. UXj^, 
UYjk) has the same vanishing entries as Hj (resp. Xj^, Yj^). 

Let us now prove that ([7j) holds for any unitary matrix. Consider U G U(A^). Choose P,D£ U(A^) 
such that D is diagonal and U = PDP~^. Set Y = P~^XP. Then [/e*^ = PDe^^P~\ The result 
follows now easily. □ 



Before we apply the last result in order to compute the differential of tr/, let us state a classical 
yet very useful lemma. 

Lemma 4.4. Let (-^fc)fcG{i,...,Af2} be a orthonormal basis ofu{N). Let A, B be elements ofM.]\f{C). 
Then the following equalities hold: 

tr{AXkMBXk) = tr(AS), (9) 

k=l 

Y^triAXkBXk) = -tv{AMB). (10) 

k=l 

Proof. 1. For A, B £ u{N), this equality multiplied by A^^ is indeed simply 

Ar2 

'^{A, Xk)u{N){B, Xk)u{N) = {AB)uiN)- 
k=l 

The general case follows thanks to the equality Mjv(C) = u(A^) © iu{N) and the fact that the 
relations are C-bilinear in {A,B). 

2. Choose l,m £ {1, . . . , A^^}. By taking A = Eji and B = Emi in the first relation, we find 

''^{Xk)ij{Xk)lm = —j;^Si^rnSj,l. 
k=l 



12 T. LEVY, M. MAIDA 

The second relation follows by developing the trace. 



□ 



Proposition 4.5. Let f : V ^ M. be a differentiahle function. Then tr/ is differentiahle and, for 
all U G U(iV) and all Y G u{N), we have 

(£y(tr/))([/) = -itr(/'([/)y). (11) 

In particular, Vf/ G U{N), ||V(tr/)(C/)||2 = ^tr(/'(C/)2). 

Proof. Since tr/ is invariant by conjugation, we have for all U,V £ U(A^) and all Y £ u(A^) the 
equality (£y (tr/))(C/) = (£yyy-i (tr/))(yC/y-^). Hence, it suffices to check ([U]) for ah Y when U 
is diagonal. In this case, the result is a direct consequence of Proposition 14.31 The second assertion 
follows from the definition of the gradient and the identity ([9]) . □ 



4.3. Lipschitz norms again. At the end of the proof of Proposition 13.11 (see Section 16. 2p . we 
will need to estimate the Lipschitz norm of a function of a unitary matrix of a special form. We 
state and prove this estimation below, although the reader might want to skip it now and jump to 
Section \5\ 

Proposition 4.6. Let f be an element of C^'^{\J). Let V,W be two elements ofU{N). Define a 
function Fy^w '■ U(A^) C by setting 

Fv,w{U) = tr {f'{UV)f'{UW)) . 

Then F is Lipschitz continuous and we have the estimate 

\\Fv,w\\up < 

Proof. We prove that Fy^w is differentiahle almost everywhere on U(A'^) and estimate the L°° norm 
of its differential. According to Proposition l4.3l we have, for all X £ u{N) and almost all U G U(A^), 
the equality 

iCxFv,w){U) =tr {V-'Df'{Lvu,Rvu){VUX)Vf'{UW)) + 

tr {f'{UV)W-'Df'{Lwu, Rwu)iWUX)W) . 

We have used the fact that i^t^of (Ue^^V) = V^^f^^^^J'{VUe^^)V. Let us focus on the first 
term of the right-hand side, the second being similar. By the Cauchy-Schwarz inequality, 

|tr {V-^Df'{Lvu,Rvu){VUX)Vf'{UW)) f < tr{M*M)tT {f'{UWrf'{UW)) , 

where we have set M = Df'{Lvu,Rvu)iVUX). 

Recah that MAr(C) is endowed with the scalar product {A,B) = NTr{A*B). We claim that 
the operator norm of the endomorphism Df'(LvU:Rvu) of M7v(C) with respect to this norm is 
bounded above by ^||/"||/,oo. Indeed, this operator is normal with respect to this scalar product, 
so that its operator norm equals its spectral radius, which is smaller than the L°° norm of D/'. 
Hence, we find 

tr(M*M)5 < |||/"||iootr(X*X)i 

It follows that 

\\^xFv,w\\l^ < 2-||/"||l°° — j^^^^ 
from which the result follows easily. □ 
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5. Convergence of the bracket 

In this section, we prove the first two assertions of Proposition 13.11 Let us first prove a funda- 
mental property of the generator of the Brownian motion on U(A^). The action of U(A^) on u(A^) 
by conjugation is an isometric action. Hence, for all V G U(A^), the processes Un and VUnV~^ 
satisfy two stochastic differential equations (see ([I])) driven by two processes in u(A^) with the same 
distribution, so that they have the same distribution. 

Lemma 5.1. Let F : U(iV) — > M 6e a Lipschitz continuous function. Let Y be an element ofu{N). 
Let t > be a real number. Then CyiPtF) = Pt{CYF). 

Proof. Since F is Lipschitz continuous, CyF is well-defined as an element of L°°(U(A^)). The result 
amounts simply to the interversion of an integration and a derivation: for all U G U(A^), 

CY{PtF){U) = E [F{Ue'^''UN{t))\ = E [F(f/[/;v(t)e^^)] 

as\s=o ds\s=Q '- 



E 



F(C/C/^(t)e^^) 

ds\s=Q 



Pt{CYF){U). 



We have used the fact that UN{t) has the same distribution as e UN{t)e . □ 



5.1. Ito formula. The following result summarizes the applications of Ito formula that we will 
use. The third assertion below implies, by polarization, the first assertion of Proposition 13.11 

Proposition 5.2. Let F : U(A^) be an integrable function. Define a real-valued martingale 

L^ indexed by [0,T] by setting, for all t G [0,r], L^{t) = E[F{UN{T))\J^N,t]- Let {Xk)k(i{i,...,N^} 
be a orthonormal basis of u{N) . Then the following equalities hold for all t G [0, T] . 

1. L^{t) = {PT-tF)iUN{t)). 

2. L^(t) = L^(0)+ / y2Cx,iPT-sF)iUN{s))d{Xk,KN)u{N)is). 

k=i 

3. {L^m= f\\{V{PT^sFmj,{s))f ds. 

Jo 

4. If F is Lipschitz continuous, then {L^){t) = / 2~][PT-s{^XkF){UN{s))]'^ ds. 

■^0 k=l 

Proof. 1. Choose t G [0,T]. Since the unitary Brownian motion has independent multiplicative 
increments, L^{t) can be rewritten as 

L^{t) = nF{UN{T))\TN^t] = nF{UN{t)UUt)UNiT))\rNA = E[F{UN{t)VN{T - t))\TN,t], 
where Vn is a Brownian motion on U(A^) with the same distribution as Un and independent of 
U]sf. The result follows. 

2. Let us apply © to the function G : [0,r] x [j{N) R defined by G{t, U) = {PT^tF){U). It 
follows from the definition of the semigroup {Pt)t>o that G satisfies the time-reversed heat equation 
^AG + dtG = 0. Hence, Ito's formula reads 

L^{t)=L^{0) + y] / {Cx,{PT^sF)){UN{s))d{Xk,KM)uiN){s)- 

3. The equality follows immediately from the equality 2 and the fact that the processes {{Xj., Kn)u(n) 
k G {1, . . . jN"^}} are independent standard real Brownian motions. 
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4. This equality follows from the previous one by applying Lemma I5.1i 



□ 



5.2. Expectation of the bracket. We can now prove the second assertion of Proposition 13. li 

RecaU that we use the notation E^^^s, U) = N'^{V{PT-s{'^^f)){U),V{PT-s{'^T^g)){U))u[N)- We will 
use the fact, which is a consequence of Jensen's inequality, that for any square-integrable function 
G : U(iV) ^ M, and for all t > 0, {PtGf < Pt{G^). 

Proof of the second assertion of Proposition \3.1[ Let / : U — > M be Lipschitz continuous. By defi- 
nition and by Lemma |5. II 

Ellf{s,U) = N'Y,iPT-s{Cx,Mmuf 

k=l 
k=l 

= iV2p^_.(||V(tr/)f )([/). 
By Proposition 14.51 and the fact that Pt~s does not increase the uniform norm, this implies that 

|i^^;^(s,c/)i<ii/'iii^. 

By polarization, the estimation of \Ej^^ {s,U)\ follows. 

Now, let us consider two independent copies Vn and Wn of the unitary Brownian motion Uj\f. 
Then, denoting by Kvj^^Wn the expectation with respect to Vat and Wn only, we have 

Ef^f{s, Uxis)) = N' Y,{Pt-s{Cx, {ir f))){UN{s)f 

k=l 

= N'Y.Ev^,wAi^x,tTf){UN{s)VM{T - s)){£x,trf){UN{s)WN{T - s))]. 

k=l 

Using successively Lemma 14.51 and Lemma 14. 4| we find 

Ef^^s, Un{s)) = Ev^,wr, [tr {f'(.UN{s)VN{T - s))/'(C/iv(s)VFjv(r - s)))] . 

Taking the expectation with respect to Un, we find finally 

E[£;^^(s, Un{s))] = E [tr {f {Un{s)Vn{T - s))f'{UN{s)WN{T - s)))] . 

Let {A, t) be a C*-probability space which carries three free mutliplicative brownian motions 
u, V, w which are mutually free. According to Theorem 1 2. 31 the family {Un{s), VN{t), WN{u))s,t,u>o, 
seen as a collection of non-commutative random variables in the non-commutative probability space 
(L°° (8) M7v(C),E <S> tr), converges in distribution to {us,vt,Wu)s,t,u>o as N tends to infinity. This 
implies in particular that for all non-commutative polynomial p in three variables and their adjoints, 
and for all s,t,u > 

E[tT p{UN{s),VN{t),WN{u))] T{p{Us,Vt,Wu)). 

N^oo 

Let us fix s G [0, T). Since ^ is a C*-algebra, there is a continuous functional calculus on normal 
elements, hence on unitary elements, and f {usVT-s)f' {usWt-s) is a well-defined element of A. On 
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the other hand, choose e > and let q{z, w) be a polynomial function in z, w and their adjoints 
such that sup^ ~ < e. Then 

|E [tr {f'{UNis)VNiT - s))f'{UNis)WNiT - s)))] - T{f {UsVT~s)f {UsWt^s))\ < 

|E [tr {f'{UN{s)VN{T - s))f'{UN{s)WN{T - s))) - tr q{UN{s)VN{T - s), Un{s)Wn{T - s))] \ 

+ |E [tr q{UN{s)VN{T - s), Un{s)Wn{T - s))] - T{q{UsVT-s,UsWT-s))\ 
+ \T{q{UsVT-s,UsWT-s)) - T^f {UsVT-s) f {UsWt~s))\ ■ 

The first and the third term are smaller than the uniform distance between q{-, •) and /'(•)/'(")) 
hence smaller than e. The middle term tends to as tends to infinity. Altogether, this proves 
that 

E[E^/{s,UNis))] T{f'{UsVT-s)f'{UsWT-s)), 

N^oo 

from which the expected result follows by polarization. □ 



6. Convergence of the variance of the bracket 
This section is devoted to the proof of the third assertion of Proposition 13.11 

6.1. A weak concentration inequality. Consider a function F : D{N) — > M. If F is Lipschitz 
continuous, then the equality ||-F||Lip = ||VF||ioo holds. The goal of this paragraph is to prove the 
following inequality. 

Proposition 6.1. Let F : D{N) ^ M. be a Lipschitz continuous function. For all T > 0, one has 

the following inequality: 

Y^t[F{Un{T))] < r||F||Lp. 

Note that this inequality is preserved by rescaling of the Riemannian metric on U(A^), that is, 
by rescaling of the scalar product on u(A^). Indeed, let A be a positive real and let us consider the 
scalar product = A(-,-)u on u(A''). Then, putting a tilda to the quantities associated with 

this new scalar product, we have on one hand d = X^d and = A~2 ||F||Lip, and on the other 

hand A = A^^A and Un{T) has the distribution of C/Af(A~^T). 

Proof. Recall the definition of the martingale L^ (see Proposition 15. 2p . The left-hand side is equal 
to E[(L^)(T)], thus, by the third assertion of Proposition 15.21 to 

r \\{V{PT^sF)){UN{s))f ds<T sup \\V{Pt^sF)\\1^=T sup \\Pt^sF\\1^. 

JO sG[0,T) se[0,T) 

On the other hand, since F is Lipschitz continuous, for all t > 0, ||-Pt-F||Lip < ||-^||Lip- The result 
follows. □ 



6.2. An estimate of a Lipschitz norm. With Proposition 16.11 in mind, we are going to study 
the Lipschitz norm of U ^ E^-^ [s, U) in order to estimate the variance of Ej^^{s, U]\[{s)). 

Proposition 6.2. Assume that f is of class C^'^(U). Then 

sup \\El^^{s,-)\\up = 0{N-'). 
se[o,T] 
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Proof. The proof rehes on the identity 

e{}Hs, Un{s)) = Kv^,wr, [tr {f\UN{s)VN{T - s))f'{UN{s)WN{T - s)))] . 

By Proposition 14.6^ the expression between the brackets is a Lipschitz continuous function of 
Un{s) for all values of Vn{T — s) and W^iT — s), with a Lipschitz norm which does not depend 
on Viy(T — s) and Wn{T — s) and is 0{N~^). Hence, the same estimate holds for the expectation. □ 



Proof of the third assertion of Proposition \3.1[ It suffices to combine Proposition 16.21 and Proposi- 
tion l6.1l to find that that sup^gjg.T) Var[£'^''^(s, Un{s))] = 0{N~^). The same result for E^^ follows 
easily. □ 



7. Other Brownian motions, on unitary and special unitary groups 

In this section, we explain how Theorem 1 2 . 6 1 can be extended to other Brownian motions on the 
unitary group and to Brownian motion on the special unitary groups. 

In this paper so far, we have considered the Brownian motian Un on U(A^) associated to the scalar 
product on u(A^) given by (X, y)„(;v) = NTt{X*Y), for any X,Y e u{N). The crucial property 
of this scalar product is its invariance under the action of U(A^) on u(A^) by conjugation. There 
is in fact a two-parameter family of such scalar products, namely a (Tr(X*y) + -^Tr(X*)Tr(y)) 
with a > and b > —1. Changing the coefficient a would simply affect the Brownian motion by a 
global rescaling of time, and we stick to a = in order to have correct asymptotics as N tends to 
infinity. However, varying b yields really different Brownian motions. We start by considering the 
limiting case 6 = — 1, which corresponds to the Brownian motion on SU(A^). 

Let us denote by su(A^) the hyperplane of u(A^) consisting of traceless matrices, which is also the 
Lie algebra of the special unitary group SU(A^) and be the linear Brownian motion on su(A^) 
corresponding to the scalar product induced by (•, •)u{n)- Let Vn be the solution of the stochastic 
differential equation 

dVN{t) = VN{t)dKl{t) - 1 (^1 - VN{t)dt. (12) 

One can check that if the initial condition is in the special unitary group, then the process Vm stays 
in it. We call it the Brownian motion on SU(iV). 

If (Yi, . . . , Y;v2_;^) an orthonormal basis of su(A^), the generator of Vn is given by 



Ar2-1 

9 2 ^ ' 



k=l 

We can modify the definition of Vjv in order to get Brownian motions with values in U(A^) : for 
all a > 0, we set 

vj^\t) = e^VN{t). 

where {Bt)t>o is a standard real Brownian motion independent of V]\f. For all a > 0, the generator 
of V^"'' is given by 



2 2 V ^InJ ' 



and we call vj^^ an a-Brownian motion on U(A^). 
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In fact, for each a > 0, the process is naturally associated with the scalar product 

{X,Y)['ll^ = NTTiX*Y) + i-^Tr(X*)TV(y) 



on u(A^). Indeed, let be the linear Brownian motion on u(A'^) corresponding to this scalar 

product. Note that it can be expressed as = Kjl^ + ^B. Then the process V^°^ satisfies the 
stochastic differential equation 

fHt) = vl-\t)dK^^\t) + ^) Vl,-\t)dt. (13) 



Equivalently, one can say that e ^V^'^ {t) satisfies the equation ([T]) in which Kiq has been replaced 

by In particular, in the case when a = 1, V^^^ has the same distribution as U]\^. 

The main feature of the Brownian motion on U ( A'^) that we have used extensively in the proof 
of Theorem 12.61 is that its generator commutes with all Lie derivatives. It is also the case for Vj\f 
and for all the processes V^"^ . 

With this definition, and following [2], one can check that for all a > 0, the process vj;^^ converges 
as tends to infinity to a free multiplicative Brownian motion. 
Let us now define a modified version of the covariance ar- 

Definition 7.1. With all the notation of Definition\2.4\ we define, for all a > 0, 



r-T 

C^T^f^d) = I T{f'{UsVT-s)g'{UsWT-s)) + ("^ " l)T{f' {UsVT-s))T{g' {UsWT-s)) ds. 





Following step by step the proof of Theorem 12.61 one finds the following result. 

Theorem 7.2. Let T > be a real number. Let n > 1 be an integer. Let /i, . . . , /n : U — > M 

be n functions of C^'^{U). Let us define a n x n real non-negative symmetric matrix by setting 

T>^\fi, . . . , fn) = {'^T^\fi, fj))i,je{i,...,n}- Then, as N tends to infinity, the following convergence 
of random vectors in M" holds in distribution: 

iv(tr/.(vt^(r))-E[tr/.(vt^(r))j), AA(0,S5?) (/!,...,/„)). (14) 

We leave the details to the reader, since every step can be adapted in a straightforward way. 
The only substantial change is in Lemma 14.41 which now will take the following form. 

Lemma 7.3. Lei (Xfc)^g|;^^ ^^2| be an orthonormal basis o/ (u(A^), (•, LetA,B be elements 

o/M7v(C). Then the following equality holds: 

triAXkMBXk) = -j^ {tT{AB) + [a" - l)tr(A)tr(5)) . (15) 

k=l 

Assume that [Xi, . . . ,Xf^2_i) form an orthonormal basis of 5u{N) endowed with the scalar product 
induced by {■,■) ^^(^j\fy Then 

7V2-1 

J2 ti{AXkMBXk) = -— {tv{AB) - tT{AMB)) . (16) 

k=l 

It is this modification which gives rise to the new covariance introduced in Definition 17.11 
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8. Behaviour of the covariance for large time 

For any fixed A^, the Markov process {Un{T))t>o converges in distribution, as T goes to infinity, 
to its invariant measure, which is the Haar measure on U(A^). In [^, P. Diaconis and S. Evans have 
estabhshed a central limit theorem for Haar distributed unitary random matrices. In this section, 
we relate our result to theirs by comparing the limit as T tends to infinity of the covariance 
with the covariance which they have found. 

8.1. Statement of the result of convergence. In order to state the result of Diaconis and 
Evans, we need to introduce some notation. 

Definition 8.1. Let H2(\J) denote the space of functions that are square-integrable on U and such 
that 



2 ._ 1 /■ \fie'n-fie''y' 



5 " 16vr^ ./[o,27r]2 sin2 ( ^ 



dipdO < oo. 



We denote by (•, •)! the inner product associated to this Hilbertian norm. 

For all / : U — > C which is square-integrable and all j G Z, we denote by aj{f) = fu f{t)e~^^^dt 

the j-th Fourier coefficient of /. One can check that / G i^a (U) if and only if '^Zj^z is 
finite and that, in this case, 

\\f\\l = Y.\^\\am'- 

2 ^ — ' 

The result of Diaconis and Evans states as follows. 

Theorem 8.2. (5.1 in [7] J For all G N, let be a N x N unitary matrix distributed according 
to the Haar measure on U(A^). Let n > 1 be an integer. For all fi, . . . , fn £ (U), /et . . . , /„) 

be the n x n real non-negative symmetric matrix defined by . . . , /„) = ( (/j, fj)i 1 . As 

N goes to infinity, the following convergence of random vectors in M" holds in distribution: 

iV(tr/,(M;v)-]E[tr/,(A%)]),g^i,,„,„} AA(0, . . . , /„)). 

In view of this result, it is natural to expect the covariance that we have introduced in Definition 
12.41 to converge, as T tends to infinity, to the covariance given by the iifa-scalar product. This is 
what the following result expresses. 

Theorem 8.3. For all n > 1 and all /i, . . . , /„ G H^{U), 

^rifi, ■ ■ ■ , fn) - — > . . . , /„). 

i — *oo 

Let us emphasize that T,T{fi, ■ ■ ■ , fn) has only been defined so far for functions in C^'^(\]). From 
this point on, we focus on extending the definition of the covariance to functions of the space 
and proving Theorem | 



8.2. The main estimate. In the sequel, {ut)t>Q, {vt)t>o and {wt)t>o will be three multiplicative 
free Brownian motions, that are mutually free. For all T > and all k G Z, let us denote by 
fJ-k{T) = t(u^) the k-th moment of ut- Recall that, since ut has the same law as n^, one has, for 
all /e G Z, the equality ^k{T) = fi-k{T). For each k > 1, according to [2], fik{T) is given by 



MT) = e-^y:i#(,M^-^ (17) 



^ l\ \l + l 

1=0 ^ 
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Lemma 8.4. For all e > 0, all T > To{e) = | log(l + |) and all k eZ, one has 

Proof. If /c = or e > ^, the inequality is trivial. Moreover, since fJ-k{T) = fi-k{T), it suffices to 
prove the inequality for k > 0. So, let us assume that e < | and A: G N*. It is easy to check that 
the expression p!7|) of /Ufc(T) is equivalent to the following: 

where we integrate over a closed path of index 1 around the origin. If we choose as our contour the 
circle of radius | centered at the origin, we get 

so that 

as expected. □ 



We will denote by Tq a real large enough such that for all T >Tq and all G Z, the inequality 
W(T)\ < e-l'^la holds. One can check that 31 is large enough but we choose Tq = 32 for reasons 
which will soon become apparent. 

For all J, A; G Z and T > 0, we define 

rT 



T [{UsVT-sy{UsWT-s) ) ds. 



(18) 



Proposition 8.5. Set Tq = 32. For all T > Tq and all {j,k) ^ (0,0), the following inequality 
holds: 



+ 



Moreover, if j / 0, then 



<^ + 2|j|roe-*(^-^o). 



lil 



(19) 



(20) 



In particular, for all {j,k) ^ (0,0), the following convergence holds : 



lim Tj,fc(T) = 5. 



1 



The proof of these estimates relies on a differential system satisfied by the functions fc. This 
differential system is a consequence of the free Ito calculus for free multiplicative Brownian motions. 
We state the form that we use, which is of interest on its own. 

Proposition 8.6. Let {ut)t>o be a free multiplicative Brownian motion on some non- commutative 
*-probability space {A,t). Let ai,...,a„ G A be random variables such that the two families 
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{ut '■ t > 0} and {ai, . . . , an} are free. Finally, choose ei, . . . , e„ G {1, —1}. Then 
^r(uf ai . . . a„) = - ^r«^ai . . . nf"a„) 

- ^ le,=ejT{ai...aj-iul')T{aj ...ai-iul') 

l<i<j<n 

+ ^ le^^e.T(ai...aj_i)r(aj...ai_i), 

l<i<j<n 

where for all 1 < i < j < n, we have used the shorthands Oj . . . Oj-i for OiU^^'^^ai+i . . . u^^'^aj-i 
and Qj . . . aj„i for aju^^^^ aj+i . . . u^"anU^^ai . . . ii^*"^aj„i. 

Proof. In [2], P. Biane showed that the free multiphcative Brownian motion {ut)t>o satisfies the free 
stochastic differential equation dut = iutdxt — \utdt, where {xt)t>o is a free additive (Hermitian) 
Brownian motion. The identity above follows from this fact by free stochastic calculus, which has 
been developed by P. Biane and R. Speicher and is exposed in [3] . □ 



ions : 



Lemma 8.7. The family iTj,k){j,k)&z2 satisfies the following system of differential equat 

Tj,kiT) = tJ-j+k{T) Tj^k{T) - - l)MTKgn{j){\j\-~l),k{'^) 

1=1 

\k\~l 

- X^d^l - H/^m(7')7-j,sgn(A:)(|fc|-m.)(r), 
m=l 

where fj^k is the derivative of the function T i-^ Tj^k{T). 

Proof. This differential system follows easily from an application of Proposition 18.61 to the expres- 
sion (llHl). □ 



Before we turn to the proof of Proposition 18.51 l^t us state some elementary properties of the 

fcT 

functions Tj^k- For all A: > 0, define the polynomial Pk by the relation fJ,kiT) = e~~ Pk{T). For 
< 0, define Pk = P-k- 

Lemma 8.8. For all j,k G the function Tj^k is real-valued and satisfies Tj^k = Tk,j = T-j-k- 
Moreover, there exists a family of polynomials {Rj,k)j,kez with rational coefficients such that the 
following equality holds : 

yj, keZ, r,- fc(T) = ^^j+k,o + e-^^i?,- fc(r). (21) 

These polynomials are characterized by the fact that for all j, G Z, Rj^ki^) = and 

lil-i 

= '^jk>oPj+k - ^ (lil - ^)PlRsgn{j)(\j\-l),k - ^ (1^1 - ^'l')PrnRj,sgn(k){\k\-m)- (22) 
1=1 m=l 

Proof. The equalities Tj^k = T-j,-k = Tk,j follow from the definition of Tj^fci using the unitarity of 
u, V, w, the traciality of r, and the fact that the families {u, v, w) and (n, w, v) have the same joint 
distribution. The fact that Tj^k is real-valued can be proved by induction using the differential sys- 
tem stated in Lemma [8. 7| or directly using the definition and the fact that (u, v, w) and {u* ,v* ,w*) 
have the same distribution. 
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The functions Rj^k defined by (j2ip are easily checked to satisfy the differential system ()22p and, 
by induction, to be polynomials. □ 

Proof of Proposition \8.5[ Since the differential equation for Tj^k expressed by Lemma 18.71 involves 
only indices {j',k') such that \j'\ + \k'\ < \j\ + \k\, we will prove the conjunction of (jl9p and (j20p 
by induction on |j| + It is understood that k = —j in ([20|) . 

The symmetry properties of Tj^k allow us to restrict ourselves to the two cases where j,k > 
and J > 0, A; < 0. We may also assume that j + k > 0. 

T 

The smallest possible value of \j\ + \k\ is 1. So, we start with Tifi(T) = T^i(T) = Te~ 2", which 

_ T T 

is smaller than e ^ for T larger than Tq. Hence, if |j| + |A;| = 1 and T >Tq, then \Tj^k{T)\ <e~^. 
This proves the result when \j\ + \k\ = 1. 

Let us consider now j and k and assume that (jl9p and ()20p have been proved for all j', /c' such 
that Ij'I + \k'\ < \j\ + |A;|. Let us first assume that j + k ^ 0. In this case, define 

Then Lemmas 18.41 and 18.71 and the induction hypothesis imply the inequality 

\Pj,k{T)\ < e^^e-— ^ + 4e^^ ^ (|j| - O^''^ ,., 

1=1 ul +1 I 

lil-i 



+ {\J\ + \k\ - l)roe^(™) Y: (IjI - l)e-'ie'^ 



+ 4e^^ 5^ (|fc| - m)e-3 " 

|fc|-i 

+ (|j| + \k\ - l)roe^(™) - l)e-^^e^^. 

m=l 

■ Tq 

Since |j| — / < |j| — ^ + |A;|, \k\ — m < \j\ + |A;| — m and e 4 < 1, we find 

lil-i 

1=1 
\k\-i 

_^4g^™T ^ g_„|:g_|j+sg„(fc)(|fc|-m)|f 
m=l 

00 

+ 2(U| + |fe| - l)2roe^(™) e-'S. (23) 

i=i 

If we are in the case where j. A; > 0, then we obtain immediately the estimate 

\P,UT)\ < { e-^^ + (8e-^^ + 2{\j\ + \k\ - ifT^e-'-^^^-^^] . (24) 

y 1 — ei2^ ^ J 

In the case where j > and A; < 0, the computation is slightly more complicated. In this case, let 
us also assume that j + /c > 0, as we have indicated that it is possible to do. Then the estimation 
of the sum over m in (|23p is the same as before, since j + sgn(fc)(|A;| — m) is positive for all values 
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of m. However, the sign of sgn(j)(|j| —l)+k now depends on /. Thus, we bound the first sum over 
I by 

j+k +00 



1=1 l=j+k+l 

iT _ . „ , T 



In the first term, we could actually have e ' 3 instead of e '12 but we are not seeking any optimality. 
In the second term, we write 



I J ~l~ I rj-i — 

and we find that the first sum over / in ([23]) is bounded by 2e — " _r ■ Finally, we have 

l-e"T3 

established that, when j > 0, A; < and j + k > 0, 

|p,.(T)| < fe-H^- + (12 e-^- + 2(U| + |^| - l)2roe-^(--o) 

\ 1 — e 12 ^ 



In view of (j24p . the last estimate holds for all values of j and k. Our choice of Tq guarantees that 
for T > To, the inequalities 

_T_ _T_ 

T e 12 e 12 1 
e"T2 + 12 < 1 and < - 

1 — e~T2 1 — e~T2 8 

hold. Hence, we find 

\P,AT)\ < (e-^^ + + 1^1 - l)^Toe-^(^-^o)) . 

Integrating the last inequality from Tq on and using the fact that ^-''^"t^^^ — liiM > lil+M^ -^yg fl^d 



|j+fc 



T 



from which it follows immediately that 



k.^mi < 4 + (Ul + |A:|)roe-^(^-^o), 

which is the expected equality. 

Let us now treat the case where k = —j. As before, we can assume that j > 0. Setting 
Pj{T) = e'-'l"^ (^Tj^^jiT) — iji^, we find, using the same estimates as before, that 

00 ^ 00 



It follows that 



so that 



|p,(r)| < 8el^-|^ V'^ +2(2U| - l)2e^™ro J^e ^12 
1=1 1=1 

\p,{T)\ < Toel^-I^o + el^-l^ + (2|j| - l)Toe^(™), 

^lUI 2; 

T 



<^ + 2blroe-^(^-^o). 



1 

which is the expected inequality. This concludes the proof. □ 
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8.3. Extension of the definition of the covariance. 

Proposition 8.9. Let f G ffz (U) be real-valued. The following properties hold. 
1. For all T>To, ^ |j/ca,(/)afc(/)r,- fc(r)| < cx). 



2. lim ^ jkajif)ak{f)Tj,k{T) 



Proof. Choose an integer n > 1. Then for all T >Tq, Proposition 18.51 implies 



J2 \3ka,{f)ak{f)T,,u{T)\< ^ ^_^|a^.(/)a,(/)|e 

\Mk\<n lil,|fc|<n 



4 



'ii+M(T-To) 



+To \mj\ + \m<^Af)ak{f)\e 

\3\,\k\<n 



\Mk\<n \im<n 

2 

IZ \ - /7-TT, . . _ / X - , .,o, , „s, _M, 



«GZ |i|,|fc|<nj+fc=i Vljl<" 

<ll/ll| (2 5:e-l'l?+2j;U|V^(^-^o) 

The first assertion follows. The second is a consequence of the second statement in Proposition 18.51 
and the theorem of dominated convergence. □ 



Proposition 18.91 above allows us to give a new definition of the covariance ctt when T is large 
enough. 

Definition 8.10. For all T > Tq and all f G H^V), we define 

'7t(/,/) = - J2 jf^'^jif)^k{f)'rj,k(T). 

j,k&Z 

Lemma 8.11. Let f be a function o/C^'^(U). For allT > Tq, the two definitions (Definition \2.4\ 
and Definition [8.10\} ofaT{f,f) coincide. 

Proof. The series X^jg^ is convergent, so that = l^|j|<„ J«j(/)e*-'* converges uni- 

formly to /' on U as n tends to infinity. Therefore, starting from Definition [2 



MfJ) = - f E jka,if)akif)i 



UsVT-syiUsWT-s)j ds. 

As the processes are unitary and X^jg^ < co, we get by dominated convergence that, for 

all T > 0, 

'7t(/,/) = - E jkaj{f)ak{f)Tj^k(T), 

j,k&Z 

as expected. □ 



24 



T. LEVY, M. MAIDA 



Theorem 18.31 is now a straightforward consequence of the polarisation of Definition 18.101 and 
Proposition 18.91 



Remark 8.12. Let us emphasize that Proposition \8.5\ implies that, for all e > and all T > Tq, 
the following series converges: 

Y^{\j\ + \k\)'-^\rjMT)\' <+oo. 

j,k&Z 



Hence, for all T > Tq, the equality 

KT{e''^,e'^) 



(i,fc)6Z2\{(0,0)} 



defines Kt as a square-integrable real-valued function on and, for all f,g ^ H^{U), one has the 
equality 

Mf,9)= I f'{e^')KTie^',e^^)g'{e^^^^ 



[0,27r] 



47r2 



We conclude this study of the covariance by showing some puzzling numerical experiments (see 
Figure [T]). It is striking on these pictures that the behaviour of the covariance (Tt(/, g) is complicated 
and interesting for small T, and much simpler for large T. It is thus not surprising that we have 
been only able to analyse ut for large T. 




Figure 1. For all fe > 1, let us define Sk{e^^) = sm{k9) and Cfc(e*^) = cos{k0). The 
pictures above are the graphs of the following functions of T for T G [0, 6]. Top left : 
(^risk, Sk) and arick, Ck) for A; G {1, ... , 8}. Bottom left : HkiT) for A; G {1, ... , 6}. 
Top center : aT{sk-,Sk+i) for k G {1,...,15}. Bottom center : aT{ck-,Ck+i) for 
k G {!,..., 15}. Top right : (Jt(s/c, Sfc+s) for k G {1,4,7,10,13}. Bottom right : 
(Tr(sfc, Sfc+2) for odd A; G {1, . . . , 13}. 
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9. Combinatorial approaches 



9.1. The differential system satisfied by the Tj^k- The differential system satisfied by the 
functions Tj^k (Lemma 18. 7p can be interpreted, at least when j and k have the same sign, in terms 
of enumeration of walks on the symmetric group, in the same vein as the computations made by 
one of us in [T7j. This is what we explain in this section. 

Fix J > 1. We consider the Cayley graph on the symmetric group 6j generated by all transposi- 
tions. The vertices of this graph are the elements of &j and two permutations cJi and o"2 are joined 
by an edge if and only if o\a^^ is a transposition. A finite sequence (ctq, . . . , cr„) of permutations 
such that (Tj and Uj+i are joined by an edge for all i G {0, . . . , n — 1} is called a path of length n. The 
distance between two permutations is the length of the shortest path that joins them. We call defect 
of a path the number of steps in the path which increase the distance to identity. Heuristically, one 
can understand the defect as follows : each time we compose a permutation with a transposition, 
either we cut a cycle into two pieces and this is a step which decreases the distance to identity, or 
we coalesce two cycles into a bigger one and this is a step which increases the distance to identity. 
The defect counts the number of steps of the second kind. 

For any o" € Sj, and two integers n, d > 0, we denote by S{a,n,d) the number of paths in the 
Cayley graph of &j starting from a, of length n and with defect d. The interested reader can find 
more details about those combinatorial objects in |17| . 

Let j,k > 1. If a £ &j and r € ©fc, we denote by o" x r the concatenation of a and r, that is 
the permutation in &j-^-k such that a x r(i) = a{i) if 1 < i < j and a x r(i) = r(i — j) + j if 
j + 1 < i < j + k. 

From Theorem 3.3 in [T7], it follows that for all T > 0, 



E 



tr([/jv(r)^>r([/jv(T)^') 



E 



ti{UN{Ty)\E tr(c/^(r)'=) 

(-T) 



n 



\J\f2d 



S{{l...j) X {l...k),n,d) 



E 



,n,d=0 



ni ,n2 ,di ,d2=0 



Inlj\l2{di+d2) 



ni\n2 



5((l...j),ni,di)5((l...fc),n2,d2) • (25) 



Moreover, for all T' > 0, we recall that all the expansions involved converge uniformly on {N, T) € 

N X [o,r']. 

Using this equality, it is for example easy to check that 



lim (E 

iV-»oo 



tr(C/;v(T)^'tr(C/;v(r)'=)] - E [tr(C/;v(ry )] E [ii{UN{Tf) 



E 

n=0 



n! 



S((l...i)x(l...J:),n,0) 



E 

ni=0 



n 



ni 



5((l...j),m,0)5((l...fc),n-ni,0) =0 



where the last equality comes from Proposition 5.3 of [T7]. Each term of the sum is indeed zero 
and heuristically, it means that a path without defect starting from (1 . . . j) x (1 . . . k) is simply 
obtained by "shuffling" two paths without defect from each of the two cycles in their respective 
symmetric group. 
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More interesting for us is the fact we can also deduce from (j25p that 

Kj^k{T) ^=^^ hm (k [tr(f/7v(Tytr(C/jv(T)'=)l - E [tr(C/iv(Ty )! E [tr(;77v(r)* 
' N^oo \ L J L J L 



oo 



e-U+kYi ^ i^5'((l . . . i) X (1 . . . fc), n, 1), (26) 

n=0 

where, a G Sj, r G and n > 1 being given, we use the notation 

S{a, ni, 1)S{t, n—ni,0)+S{a, ni, 0)5(r, n — ni, 1) 



S'{aXT,n,l) = S{aXT,n,l)-Y^ (^^^ 



ni=0 

Thus defined, S'{a x r, n, 1) is the number of paths of length n starting from a x t such that the 
unique step which increases the distance to the identity is the multiplication by a transposition 
which exchanges an element of {1, ■ ■ ■ ,j} with an element of {j + 1, . . . , j + k}. Thus, heuristically, 
the unique step which is a coalescence is a coalescence between a and r. 
Our goal is now to show the following combinatorial identity 

Proposition 9.1. For any integers j, k > 1, and n > 0, we have 
S'{{l...j)x{l...k),n + l,l) = jk S{{l...j + k),n,0) 



fc— 1 n / \ 

■ ■ ■ m),q,0)S'i{l ...j)x{l...k-m),n-q,l). 

m—^ n—n V^/ 



m=l 9=0 

The combinatorial interpretation of this identity is the following : let us consider a path of 
length n + 1 from (1 . . . j) x [1 . . .k) whose unique step increasing the distance to identity is a 
true coalescence between the two cycles. The first step of such a path can be of three kinds, 
corresponding respectively to the three terms of the right hand-side : 

- either it coalesces the cycles, creating a (j + A;)-cycle, and this can be done by choosing an 
element in each cycle. Then the path can be completed by any path of length n without 
defect from a (j + /e) -cycle. 

- either it cuts the cycle (1 ■ ■ ■ j) into two cycles, one of length I that will then be cut p times 
without being affected by the coalescence and another of length j — I which contains the 
element which will be exchanged with an element of {j + 1, . . . , j + fc} during the coalescing 
step. 

- either, symmetrically, it cuts the cycle [1 . . .k). 

We will hereafter propose a rigorous proof of this identity through the free stochastic calculus 
tools introduced above in the paper. 

Proof. Let the integers j,k > 1 and the real T > be fixed. If we consider the quantities K,j^k{T) 
as defined in (|26p . if we denote, for any r G Z, by /r- : U ^ C the function given by /r(-z) = z^, 
then, from Definition 12.41 and Theorem 12. 6^ we get Kj^kiT) = o'T{fj,fk) and from (fT8|) . it can be 
reexpressed as Kj^kiT) = —jk Tj^kiT). Now, from Lemma \87l\ we get immediately 

• + k ''^'^ 
Kj,k{T) = -jki2j+k{T) — K,j,k{T) - j'^^ii{T)aj_i^k{T) - A; E f^m{T)Kj^k-m{T), 

1=1 m=l 
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SO that we get immediately the anounced result, as we know from p!7j that, for any r G 

//,(r) = e-^^5^^^5((l...r),n,0) 

and from (1261) that 



n=0 



^j,kiT) = -^-^n^,k{T) - f; ^^^'((1 . . -i) X (1 . . . fc), n + 1, 1). 

n=0 



□ 



9.2. Mixed moments of special unitary matrices. In principle, any computation involving 
functions invariant by conjugation on the unitary group can be performed by using harmonic 
analysis, that is, the representation theory of the unitary group. In this section, we use this 
approach to prove the following formula, which yields for each > 3 an explicit expression for the 
covariance of traces of powers of the Brownian motion on SU(A^). With the help of Section [3 it is 
easy to deduce the analogous result for the Brownian motion on U(A^). 

Proposition 9.2. Let N > 3 be an integer. Consider, on SU(A^), the Brownian motion {VN{t))t>o 
associated with the scalar product (-'^i ^)su(Af) = NTi:(X*Y) on 5u{N). Let n and m be positive 
integers. Assume that N > n + m + 1. Then 



E 



"I , r , \ t n(n — l)-'rmhn — l) t 

T¥(V^(t)")Tr(y^(t)-) = n6n,m + (-l)"+-e-("+'")5 n 2" 



1 1 

EE 

ri=0 r2=0 



ri J\ n J \ r2 J\ m 



{N + ri+r2 + l){N -n-m + ri+r2 + l) 
{N -n + ri+r2 + l){N -m + ri+r2 + l) 

The basic strategy for the proof is to expand the heat kernel and the traces in the basis of Schur 
functions, and then to use the multiplication rules for Schur functions and their orthogonality 
properties. The multiplication rules are expressed by the Littlewood-Richardson formula and they 
are rather complicated. Fortunately, in the present situation, the Young diagrams which occur are 
simple enough for the computation to be tractable. 

Let us recall the fundamental facts about Schur functions. Details can be found in [8]. A Young 

diagram is a non-increasing sequence of non-negative integers. If A = (Ai > A^ > 0) is such 

a sequence, we call k the length of A and denote it by ^(A). The set of Young diagrams of length at 
most k is denoted by . We draw Young diagrams downwards in rows, according to the convention 
illustrated by the left part of Figure [21 

The Schur function s\ is a symmetric function which, when evaluated on strictly less than £(A) 
variables, yields 0. Whenever ^(A) < A^, the function s\ is well defined and non-zero on SU(A^). Its 
value s\{Im) at the identity matrix in particular is a positive integer, which is the dimension of the 
irreducible representation of SU(A^) of which s\ is the character. Another number attached to A will 
play an important role for us, which is the non-negative real number c(A) such that /S.s\ = —c{\)s\. 

It happens that distinct Young diagrams yield the same function on SU(A^) : if A and /i are 
Young diagrams such that i{X),£{iJ,) < N, then sx = if and only if there exists I £ Z such that 
X = fi -\- {I, . . . , I) = fi -\- . In fact, if px and are the representations of U(A^) corresponding to 
A and fi, then px = det*^^ and the restrictions of these representations to SU(A^) are equal. 
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n — r 



t U 



Figure 2. The Young diagram on the left is (7, 6, 3, 3, 1), which we also denote by 
76 3^ 1. The diagram on the right is rjn,r = (n — 1) l*". 



Finally, we need to use the decomposition of the heat kernel and the function U i— > Tr(C/") in 
terms of Schur functions. For the latter, we introduce a class of Young diagrams called hooks. For 
all n > 1 and all r G {0, . . . , n — 1}, we define 

r]n,r = (n - r, 1, . . . , 1) = (n - r) V, 

r 

which is depicted on the right part of Figure [2l 

The heat kernel at time t on SU(A^) is the density, denoted by Qt : SU(A'^) M, of the distribution 
of VN{t) with respect to the Haar measure. 

Proposition 9.3. Choose N > 1 and U G SU(A^). Then the following equalities hold. 

n.-l 

1. For all n>l, Tr([/") = J^(-l)"s^„,,(f/). 

r=0 

2. For all t > 0, Qt{U) = ^ e-'-^^sx{lN)sx{U). 

The proof of the first equality can be found in |19j , the proof of the second in [18] . The expectation 
that we want to compute in order to prove Proposition 19.21 is thus equal to 

n— 1 m— 1 „ 

The multiplication of Schur functions is governed by the Littlewood-Richardson formula, which 
describes a non-negative integer A'^^ « foi^ each triple of Young diagrams a, f3, 7, in such a way that 

7 

Using these coefficients, the integral above can be rewritten as 

JsUiN) 7 JSU(Af) 

= Yl ^lvn,r^ Yl l7='?m,r2+«^ 

7 />0 
l>0 
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Thus, we need to compute 

n — 1 m— 1 

j_]N 

ri=0r2=0 />0 

It turns out that a shghtly more general computation is simpler to perform : we compute the 
Littlewood-Richardson coefficient Na^-q^.r foi^ oi) ^i^d all n, r. Let us introduce some notation. 

Let a = (ai, . . .) and (5 = . . .) be two Young diagrams. Set \a\ = ai and = Yli l^i- 
We assume that a C /3, that is, ai < (3i for all i. Then we denote by /3/a the set of boxes of 
the graphical representation of (5 which are not contained in a. We say that a subset of f3/a is 
connected if one can go from any box to any other inside this subset by a path which jumps from 
a box to another only when they share an edge. 

We denote by k{P/a) the number of connected components of 13 /a. Also, we define v{f3/a) as the 
number of boxes of P/a which are such that the box located immediately above also belongs to P/a. 
Alternatively, this is the number of distinct occurrences of the motif formed by two consecutive 
boxes one above the other in P/a. 

Our main combinatorial result is the following. 

Proposition 9.4. Let a and /? be two Young diagrams. Let r]n,r be a hook. Then Na,rin r non- 
zero if and only if the following conditions are satisfied : a C P, \P\ = \a\ + n, P/a contains no 
2x2 square, and v{P/a) < r < v{P/a) + k{P/a) — 1. In this case, Na,r]„,r = it^-vlji/a)) ' 

Proof. According to the Littlewood-Richardson rule, Na,n„ r is the number of strict expansions of 
a by rjn,r which yield P, that is, the number of fillings of P/a with the boxes of r]n,r such that the 
following conditions are satisfied: 

1. for all s > 1, the union of a and the boxes of P/a filled by the first s rows of rjn^r is a Young 
diagram, 

2. no two boxes of the first row of ?7„,r are put in the same column of P/a, 

3. if one goes through the boxes of P/a from right to left and from top to bottom, writing for 
each box the number of the row of r]n,r from which is issued the box which has been used to fill 
it, one obtains a sequence which starts with 1, and in which all other numbers 2, . . . , r appear, not 
necessarily consecutively, in this order. 

It is important to notice that, according to the third rule, a strict expansion of a by a hook 
which yields P is completely characterised by the set of boxes of P/a which are filled by boxes 
issued from the first row of the hook. We say for short that these boxes of P/a are filled by the 
first row. 

The first two conditions a C P and \P\ = \a\+n are obviously implied by this rule. A less trivial 
implication is that there cannot exist a strict expansion if P/a contains a 2 x 2 square. Indeed, 
by the first two rules, the bottom-left box of the square cannot be filled by the first row and the 
bottom-right box must then be filled with a boxed issued from a strictly lower (in the graphical 
representation) row of rjn^r- This contradicts the third rule. 

Let us assume that P/a contains no 2 x 2 square. Then each connected component of P/a is a 
"snake" (see Figure [3]). 

Any box of such a snake which has a box on its right must be filled by the first row. These boxes 
are the white boxes in Figure [3l Any box located below a white box cannot be filled by the first 
row. These boxes are the grey boxes in Figure [3l Only one box is not in one of these two cases, 
the top-right box of the snake. In the topmost connected component of P/a the third rule implies 
that this box must be filled by the first row. 

Finally, if the first three conditions are satisfied, then P/a contains one box in each connected 
component, except the topmost one, which can either be filled by the first row or not. The minimal 



30 



T. LEVY, M. MAIDA 




Figure 3. White boxes must be filled by boxes issued from the first row of 77n,r- 
Grey boxes cannot. The black box may or may not, except if this snake is the 
topmost connected component of /3/a, in which case it must also be filled by a box 
issued from the first row of ?7n,r- 

number of boxes which are not filled by the first row is the number of grey boxes, which we have 
denoted by v{P/a). This is the minimal value of r for which there exists a strict expansion of a by 
T]n,r which yields /?. Moreover, for this value of r, the expansion is unique, since the boxes filled by 
the first row are completely determined. Similarly, the maximal value of r is v{f3/a) + k{(5/a) — 1. 
For r between these two bounds, there are exactly as many expansions as there are choices of which 
snakes have their top-right box filled by the first row. There are thus (^^^(^/q)) such expansions. 
□ 



Corollary 9.5. Let a and (5 he two Young diagrams. Choose n>l. Then 

r=0 

if a d (3, \(3\ = \a\ +n, j3/a contains no 2x2 square and is connected. Otherwise, it is equal to 0. 

Proof. If the first three conditions are not satisfied, then Na,ri„ r — ^ fo'^ all r = ... n — 1. Let us 
assume that they are satisfied. Then, by the previous proposition, the sum above is equal to 

v{l3/a)+k{l3/a)-l 



\r-v{(3/a))' 



r=v{j3ja) 

which is equal to unless k{(3/a) = 1. In this case, only one term of the sum is non-zero, for 
r = v{[3/a). □ 

We apply now this result when (3 is of the sum of a hook and a rectangle. 

Lemma 9.6. Consider n > m > 1, r2 £ {0, . . . , m— 1}, and N > m + n. For all ri £ {0, . . . , n— 1}, 
define 

<r2,n,ri = (n - n + m - Ta) (n - n + ly^ (n - nf-^^-^-^ (n - n - 
Then, for all A G Nf^'^ and all I > 1, 



,r,„,.,+/- ^ / (-1)""' z/ / G {1, . . . , n} and A = A^_,^_„_„_,, 
otherwise. 

ri=0 ^ 

Moreover, when I £ {!,..., n}, the only non-zero term of the sum is the term corresponding to 
ri = n — I. 

Finally, if n = m, then ^x^l^l = ^ if = f2 o.nd A is the empty diagram, and otherwise. 
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N 



771 - 7-2 



r 



ri 



n — Ti 



Figure 4. The diagram X^^r2,n,rv 



Proof. Let us first consider the case n > m. In this case, according to Corollary l9.51 in order for the 
sum to be non-zero, A must be a Young diagram of length at most — 1, contained in 7?m,r-2 + 
such that {rim,r2 +^^)/A contains no 2 x 2 square and is connected. Since n > m, I must be positive, 
so that the diagram r/m,r2 + has length N whereas A has length at most — 1. Thus, the A^-th 
row of {f]m.r2 ~l~ ^^)/^ Tiot empty, it has actually length /. In particular, \f]m.r2 + - |A| > I. 
If I > n, all the Littlewood-Richardson coefficients appearing in the sum are zero. Otherwise, if 
I < n, there is exactly one way to choose A a subdiagram of r/m.r2 + such that all conditions are 
satisfied : it is A = X^^^^^^^^i-i- 

When n = m, nothing changes for / > 1. However, the sum may be non-zero even for / = 0. The 



diagram A must be the empty diagram and it is easy to check that A^j 



i^ri ,T2 ■ 



□ 



We can now go on to compute (j27p . We find the following result. 

Proposition 9.7. Let N , n and m be three positive integers. Assume that n > m and N > n+m+1. 
Then 



E 



TTiVjvitrmvNity 



n. 



71—1 771—1 

ri=0 r2=0 



(In)- 

m.rn .n.r-- 



Proof. We have 

71—1 771—1 711— 1 71 

_i_jN 



ri=0 r2=0 



r2=0 1=1 
71—1 m— 1 

ri=0 r2=0 



\_\JV 

m,r2 .n,n — l 



ri+r2 I 



\_\JV 

m,r2,Ji,7"i 



The claimed equality follows easily. 



□ 



In order to prove Proposition [921 there remains to compute c(A^ 



This 



^m,r2,n,n) and Saa;,,2,„,,^ (In) 
is by no means complicated but slightly tedious. We recall the general formulae, give the results in 
this particular case and invite the reader to check them if s/he feels inclined to do so. 
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Lemma 9.8. Consider n,m > 1, ri G {0, . . . , n — 1} and r2 € {0, . . . ,m — 1}. Then the following 
identities hold. 



n{n — 2ri — 1) m{m — 2r2 — 1) {n — m)^ 

m,r2,n,ri> = H \- m -\ — . 



{N — ri — r2 — 1){N + n + m — ri — r2 — 1) 
- " ^ ' (jv + n — ri — r2 — 1)(-/V + m — ri — r2 — 1) 



m,r2 ,Ti,r-|^ 

n — 1\ f N + n — ri — 1\ [ m — 1\ f N + m — r2 — 1 



ri J \ n J \ T2 J \ ^ 

N 



Proof. The general formulae are the following : for all a G , one has 

\'i=l l<i<j<N J \i=l 



OLi 



on one hand and, using the notation A(A) = ni<j<i<Af('^« ~ ^i) ^ ~ ~ ~ 2, . . . , 1, 0), 

A(a + 5) 
A(a) 

on the other hand. □ 



Proposition 19.21 now follows easily from Proposition 19.71 and Lemma 19.81 
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